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AN IMPROVED TECHNIQUE FOR DETERMINING REFLECTION 
FROM SEMI -INFINITE ATMOSPHERES WITH LINEARLY 
ANISOTROPIC PHASE FUNCTIONS 

Clifford L, Fricke 
Langley Research Center 

SUMMARY 

A solution to the problem of reflection from a semi -infinite atmosphere is presented, 
based upon Chandrasekhar's H-function method for linearly anisotropic phase functions. 

A modification to the Gauss quadrature formula which gives about the same accuracy with 
10 points as the conventional Gauss quadrature does with 100 points is developed. A 
computer program achieving this solution is described and results are presented for sev- 
eral illustrative cases. 


INTRODUCTION 

In 1949 Chandrasekhar (ref. 1) published solutions to the radiative transfer problem 
for certain special phase functions. The solutions involved his H-functions. These func- 
tions depend on several parameters and cannot be simply defined, but can be computed by 
using methods described in reference 1, These procedures were quite cumbersome with 
the computing methods available in that day, especially when it is realized that the Gauss 
quadrature requires so many points. A modification of the Gauss quadrature introduced 
herein achieves better accuracy than that obtained by Chandrasekhar and with far fewer 
points. This is an important result since it can then be applied to more complex prob- 
lems where the number of Gaussian points required for acceptable accuracy might exceed 
the capacity of even the most modern computers. 

The computer program presented obtains solutions using the linearly anisotropic 
phase function Wo(l + x cos 0) for the problem of reflection in any direction from a 
semi-infinite atmosphere for radiation incident in any direction. The modified or 
composite -Gauss quadrature is used and achieves the same accuracy with about one- 
tenth the number of points as the usual Gauss quadrature. 



SYMBOLS 


Algebraic 

FORTRAN 


ai 

A(I) 

weight at discrete values of Mi 

c 

C 

quantity which is a function of moments of H-function, 
defined in equation (10) 


DT 

increment given to XI 

E 

E 

exponent for division of intervals 

F 

F 

1/tt times net flux of incident radiation or 1 /tt times 
rate of flow of radiant energy across a surface element 
per unit area 

H(fXi) 

H(I) 

H-function ,at discrete values of pj 

H(0)(|Lti) 

HMO(I) 

zeroth order H-function at discrete values of 

H(D(Mi) 

HMl(I) 

first order H-function at discrete values of p^ 

I 

ET 

intensity of radiation in a certain direction or radiant 
energy per unit time transported across an element of 
area normal to radiation in an element of solid angle 


EO 

that part of intensity independent of <p 


El 

that part of intensity dependent on (p 


ETXO 

ratio of reflected to incident radiation or I/Pq 

k 


independent variable of S^^^(k), defined in equation (6) 

ko; 

XKO(I) 

characteristic roots for zeroth order ^ 


XKl(I) 

characteristic roots for first order 
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Algebraic 


FORTRAN 


i 

L 

order of term in Legendre pol 3 momial expansion 

m 

M 

number of subintervals and equal to n/2 

n 

N 

number of discrete directions (positive p only) in quad- 
rature or integration sum 

P 


phase function 

R 

R 

total reflectance 

S 

S 

a function of k which gives characteristic equation when 
set equal to zero 

SN 

SN 

quantity showing sign of S on left side of a pole 


SS 

derivative of S with respect to k 

X 

X 

anisotropic parameter 


X(J) 

boundaries of subintervals 

«o 

ALPHO 

zeroth moment of -function 

«1 

APLHl 

t 

first moment of h(®) - function 

T 


optical depth 

Q 

THET 

angle relative to vertical 

0 


angle between incident and scattered ray for a single 
scattering 


PHI 

azimuthal angle 


PSI 

characteristic function 
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Algebraic 


FORTRAN 



XI 

= cos 9 

Mi 

XM(I) 

discrete value of p 

Mo 

XIO 

direction of incident radiation 

U>0 

WO 

albedo 

Subscripts: 



i 

I,M ■ 

= 1, 2, . . ., n 

j 

J 

= 1, 2, . . ., m 

0 

0 

= direction of incident beam 


THE EQUATION OF RADIATIVE TRANSFER 

Figure 1 illustrates a beam of light Iq incident on a plane atmospheric surface at 
a polar angle of 9q. In general, the intensity of the scattered radiation is a function of 
T the optical depth in the negative z direction, p the cosine of the polar angle 9, 
and 4> the azimuth angle with respect to the X-axis, as illustrated in figure 1 which 

z 



Figure 1.- Illustration of incident beam Iq and a scattered ray I(t,h,0). 
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shows the ray Rays in the downward direction have negative jx, so that the 

incident beam has a direction cosine -Mo where Mo is a positive number. The equa- 
tion of transport appropriate for the diffuse reflection and transmission of radiation 
(ref. 1, p. 22) is given by 

^ i(j fj, (j,) - i p(m,</>; M’,<^') I(t,M’,</>') dM' d^’ 

dr 4n Jg 

- 1^ e-'^/Mo -Mo,<^o) 

The term on the left is associated with the change of intensity as the depth is pene- 
trated. The first term on the right is associated with the absorption. The second term 
is the contribution from light scattered in the direction (m, 0) from all directions (m’,^’). 
The last term can be identified with the incident beam which acts as a forcing function 
and drives the diffuse radiation. Solutions are required which satisfy the following 
boundary conditions at r = 0 and t = <» (semi -infinite atmosphere): 

I(0,-M,</>) = 0 (0 ^ M S 1) (2) 

lim I(t,±m,0) = 0 (3) 


The phase function of present interest is the linearly anisotropic expression 

p(cos 0) = o)o(l + X cos 0) (4) 


where 


cos 0 = mm’ + \/l - M^ \/l - M'^ cos i<p' - (p) 


When X = 0, the isotropic case is obtained. Other phase functions have been considered 
by Chandrasekhar but will not be considered herein. 

For phase functions which can be expanded as a series of Legendre polynomials, 
the solution can be expressed in the form 


I(t,M,<^) =^I^^\t,m) cos £(0 - < Po ) 
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where It corresponds to the order of the Legendre polynomial and has the values 0 and 
1 for the phase function of present interest. 

THE METHOD OF SOLUTION 

Chandrasekhar's method of solution assumes the radiant field to be divided into 
2n streams in the directions Mi where i = ±1, ±2, . . ,, ±n and is often referred 
to as the discrete -ordinate method. The integral equation is then replaced by a system 
of 2n linear equations. Associated with each direction Mi is a weight aj appropri- 
ate for an integration formula based on the division of the interval (-1,+1). For example, 
in this approximation, the integral 



can be expressed as the summation 


±n 

I 

i=±l 


I(Mi) 


By increasing the number of directions, one would expect to approach the exact 
solution as a limit, A division according to the Gauss method is used by Chandrasekhar; 
however, it will be shown that another division is better suited to this particular problem. 

Chandrasekhar (ref. 1, p. 127) expresses the solution to the radiation equation in 
terms of H -functions indicated by the following equation: 


h^^)(m) = 



K=1 (1 + ko;M) 


(5) 


where i refers to the order of the H-function (which corresponds to the order of the 
Lengendre polynomial expansion of the phase function) and k^, are the roots of the 
characteristic equation (when S = 0) given by 


s(^)(k) = 


1 b - 2,,.2 


( 6 ) 
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The characteristic function depends on the particular phase function. For the 

phase function in equation (4) 


^f(O) = 1 + x(l - 


(7) 


and 


^(1) = 1 xo)o[l - fx2] (8) 

Note that when x = 0, this case reduces to isotropic scattering and is not 

needed. 

REFLECTION IN TERMS OF THE H -FUNCTIONS 

The solution to the equation for the phase function a>o(l + x cos 0) in the case of 
the semi-infinite atmosphere (ref. 1, p. 138) is given by 

I(O,M,<^;fio0o) = H(0)(/i.o)[l - c(fi + iiq) - x(l - Wo)mMo] 


(9) 


+ X 


\/l - \/l - h(1)(m) h(1)(mo) cos (<f>o - 


where 


C = XO)o(l 


- Wo) 


^1 

2 - woQfo 


0!o = 



and . 



( 10 ) 


( 11 ) 


( 12 ) 


7 



Another quantity of interest is the total reflectance in the vertical direction which 
is defined by 


R(‘^o,x,Mo) 


' - nl p27T 

— i— \ \ I(O,m,0)m dM d0 

^loffF v)q Jq 


(13) 


Thus, 


R(coo,x,Mo) = ^ h(0)(mo) h(0)(m) 


1 - x(l - a)o)jijio 
M + Mo 


- c 


IX dll 


(14) 


This integration and those for the moments will be computed by using the summa- 
tion at the same discrete values of jx and the same weighting as’used in the calculation 
of the H-functions (eq. 5). 

PRELIMINARY DISCUSSION OF COMPUTATIONAL METHODS 

The main difficulties encountered in Chandrasekhar's H -function method are due to 
the large number of discrete ordinates required to give acceptable accuracy and the cor- 
responding large number of roots of S obtained in equation (6), These difficulties are 
greatly diminished by a modification of the Gauss quadrature. Discussions of this 
modification and properties of the S-equation used in the computer program to find the 
roots follow. 


Integration Intervals 

With modern computers it is possible to calculate to much higher orders of dis- 
crete ordinates than did Chandrasekhar who used an iterative scheme and interpolation to 
calculate the higher orders (ref. 1, p. 123). The computer program presented herein was 
used in order to test the variation of the H-f unction for the Gaussian integration as a 
function of the number of points up to 100 for isotropic scattering (i.e., oJq = 1) corre- 
sponding to table XI, p. 125, of reference 1. Figure 2 shows plots of h(®) for 
IX = 0.05, 0.5, and 1.0 as a function of the number of Gaussian points. The slowest con- 
vergence occurs for p = 0.05 where it is noted that even with 100 points an error of 
about 0.0001 occurs. 

The Gauss method will exactly integrate polynomials of degree less than 2n where 
n is the number of Gaussian points. The difficulty is that the polynomial expansion of 
the integrand over the entire interval is of a higher order than that which can be practi- 
cally used. Thus, it seems feasible to obtain a better approximation by breaking up the 
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(c) n = 1.0. 
Figure 2.- Concluded. 


entire interval into several subintervals, each of which can conceivably be expressed 
very accurately by relatively low -order polynomials. This procedure has been called a 
composite quadrature (ref, 2). The Gauss method places more points near the /i = 1 
end of the integration interval, but it will be shown that better accuracy is obtained if 
more points are placed near the /i = 0 end. A scheme which gives excellent accuracy 
with few points is to break the interval (0,1) into a number of subintervals placing more 
points where needed and then within each subinterval use a two-point Gauss quadrature. 
A convenient way of doing this is to space the intervals according to the formula 


Mj = 



(j = 1, 2, . ., m) (15) 


where fij is a boundary between subintervals and m = n/2. 

The weight for each point in the subinterval is equal to half the width of that sub- 
interval (ref. 3, pp. 887 and 916): 


Mj 

^2j = a2j-l= — 


-Mj.i 

2 
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The corresponding two discrete values of (x in the subinterval between pj and 
are given by the following equations (ref. 3, p. 916); 


^2j-l 


Mj + Mj-l fij - fij-l 
2 ' 2\/3 


and 


M2j = 


fJ-j + Mj_i Mj - Mj_l 

2 2y/3 


Equal subintervals correspond to E = 1 and would be expected to be superior to 
a Simpson Rule. For E=2 the intervals are bunched at the ii = 0 end. This scheme 
is used in subroutine SNTRVL presented later. . 

The program was run for various values of the number of points N. Fig- 
ure 2 shows the convergence of as a function of N for the cases E = 1.0, 1.5, 

2.0, and 2.5. An isolated point for E = 0.55 is also shown. With E = 2.0 and 
N = 10, the largest error is about 0.0001 for the worst case (p = 0.05) and is comparable 
with the 100-point Gauss method. With 20 terms the error is less than 1 x lO"^. It may 
be noted that the error in h( 0) (0.05) in reference 1 is 0.000225. With 100 terms and 
E = 2, the following values are obtained with eight decimal accuracy: 


h(0)(0.05) = 1.136 574 83 


h(0)(0.5) = 2.012 778 77 


h(0)(1.0) = 2.907 810 53 


Roots of the S- Equation 

Since equation (6) for S must be solved for the n values of k which give 
S = 0, it is instructive to examine the properties of S as a function of k. It may be 
noted that 4^^®) > 0 for all values of x, p, and Wq, but always has the same 

sign as x. Thus, 




^=0 

O 

II 

II 

o 

lim s(^)(k) = lim 
k-0 k-0 

-l + 27T^aj^W(p.) 

1 

o 

II 

■ V 

o 



.<0 

(£ = 1) 


(16) 
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and 


, . lim ,s(^)(k) = ±00 
‘ k-l/Mi 


Thus, there are n poles located at kj = In addition, 


(17) 


lim S = -1 (18) 

k— 00 


(_x)^ 

Consider also the quantity SN = 
the sign of and indicates the sign of S 


where |x| > 0; this quantity always gives 
just to the left of a pole. 


^ These parameters are shown schematically in figures 3(a) and (b). These proper- 
ties will' be used in the computer search for roots. 



(a) Locations when SN = +1 and either Z = 0 for any x or £ = 1 for x > 0, 



1/y^ 1/Pn-i I/M3 1/^2 I/VI3 


. . (Td) Locations when SN = -1 and ^. = 1 for x < 0. 

Figure 3.- Schematic locations of roots and poles of S as a function of k. 
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DESCRIPTION OF COMPUTER PROGRAM 


The main program CHSKR calls subroutine INTRVL and subroutine HFUNC (which 
in turn calls subroutine ROOT), The programs and subprograms are listed in appendix A, 
Flow charts of the main program and subprograms are shown in appendix B. The sub- 
routines will be discussed first. 

Subroutine INTRVL (N,XM, A) 

This subroutine computes the n discrete values of juj storing them in XM and 
the weights aj storing them in A(I). 

Subroutine ROOT(Xl,DT,SN,N, XM, A, WO, L,X) . ' • v , 

This subroutine finds a root of the* S-fvmctipn which. as known, to lie4n the Interval 
between the two points XI and X1-2*DT. The first trial root is . XI .incremented by 
DT (where DT is equal to half the interval). The values of S and SS (the deriva- 
tive of S) are computed first; then DT is halved. If the new value of S has changed 
sign, then the sign of DT is changed. 

By Newton's method, the next increment would be S/SS. However, if the Newton 
increment is smaller than DT, then DT is set equal to S/SS. This prevents having 
an increment which is so large that the next trial value is completely out of the interval; 
such an increment can occur in the Newton method if the trial point is too far from a root. 

When the increment DT is small compared with the value of the root XI, then 
the search is over and the root XI , is returned. 

The ejq)ression for PSI in the subroutine gives the correct expression for 
whether L is 0 or 1, 


Subroutine HFUNC(WO,X,K,N,XM,A,XK,H) 

This subroutine computes the H -function for either L = 0 or 1 and stores the 
value in H(M) corresponding to each jii or XM(M) value. The roots ki of the 
function S are stored in XK(M). 

The first statement sends the computer to statement 30 if X = 0 and L = 1, 
because then h( 1) does not exist. In order to indicate this condition in a convenient 
manner, DO-loop 35 sets each root kj equal to l/pj so that subsequent calculation 
will always give H = 1. (The actual value given to H is unimportant because x = 0 
in eq. (9) will always make that term zero.) 

The sign of SN indicates the sign of S just to the left of a pole. Each root 
except the last one (when J = N) lies between the poles l./XM(J) or XI and 
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1/XM(J+1) so the increment DT is set equal to (Xl-l./5^(J+l))/2. Subroutine 
ROOT then finds this root and it is then stored in XK(J). 

The last root (for J = N) can be equal to 0, less than 1/XM(M), or greater than 
1/^(1). If SO, the value of S at k » 0, is greater than or equal to 0, the root 
must be located at 0 (slight numerical errors may cause SO to be slightly positive). 
However, if SO differs in sign from SN, then the root lies to the right of 1/XM(1), but 
will always be less than 2/XM(l), Otherwise, the root lies between 0 and 1/XM(N). 

Program CHSKR 

The main program calculates all quantities for x = 1, 0, and -1, and for ' Wq = 1.0 
and 0.8. 

Subroutine HFUNC is called first with i = 0 and then with jf = 1 in order to 
compute the zeroth and first-order H-functions. Roots for the zerpth order character- 
istic equation are stored in XKb(I), while the H-function at discrete values of is 
stored in HM 0(1); the first-order quantities are stored in XKl(I) and HMl(I). 

The variable Mo is XIO, The quantities h(0)(mo) and h(I)^Mo) are called 
HOO and HIO in the program. 

The reflected intensities ET are broken into two parts: EO is the part of the 
intensity which is dependent on 9 and El is the part which is dependent on <p. 

ILLUSTRATIVE EXAMPLES 

Appendix C shows the results of running the prc^ram with N = 30 for the cases 
X = 1.0, 0.0, and -1.0; Wq = 1.0 and 0,8; and Po = l-O, 0-8, and 0.6. The total running 
time on a Control Data 6400 computer system was about 21 seconds. Other cases can 
easily be obtained by slight changes in the program. The h( 0) -functions for the cases 
X = 0.0 with ojo = 1.0 and 0,8 may be compared directly with reference 1, table XI, 
page 125; the case x = 1,0 with wq = 0.8, may be compared with reference 1, table XVI, 
page 139. The -function for the cases xo)© = 1.0, 0.8, -0.8, and -1,0 may be com- 
pared with reference 1, table XVm, page 141. These comparisons are shown in tables I 
and n for selected values of fi. The agreement is excellent in all cases, except for 
small values of fi where a difference of 1 or 2 exists in the fourth decimal. 

The reflection is shown for angles of 9 from 0° to 180° in steps of 10° and for 
angles of 0 = 45°, 90°, 135°, and 180° where appropriate. The reflection for the cases 
X = +1.0 and -1.0 with Wq = 1.0 and fiQ = 0.8 agrees with reference 1, figure 11, 
page 148. The case for x = 1.0, Wq = 0,8, and Mo = 0.6 agrees with reference 1, fig- 
ure 12, page 149. 
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CONCLUDING REMARKS 


A computer prc^ram has been presented which computes the reflection from a 
semi -infinite atmosphere for the linearly anisotropic phase function. This program uses 
the method of Chandrasekhar, but with a new modification to the conventional Gauss quad- 
rature formula. This departure from the traditional quadrature reduces the number of 
points so that answers for any set of parameters can be easily obtained with high accu- 
racy and very short computer running time. The computer program presented can be 
easily modified for other phase functions which are capable of being solved by this 
method. This composite Gauss quadrature should also be valuable in more complex 
problems where the number of points must be limited. 

Lai^ley Research Center 
National Aeronautics and Space Administration 
Hampton, Va. 236615 
July 14, 1975 



APPENDIX A 


LISTING OF COMPUTER PROGRAM 

This appendix contains a computer program for computii^ the reflection from a 
semi-infinite atmosphere for the linearly anisotropic phase function. This program uses 
the method of Chandrasekhar, but with a modification to the conventional Gauss quadra- 
ture formula. 


PRUoRAM CHSRK( OUTPUT) 

♦.♦♦♦♦♦this program COMPUTfcS THL Rt^LECTlU!M ET iN ANY DIRECTION (THETA, PHI) 

» ERUiM A SEiMl-lNE IMTc ATMOSPHERE rlAVING A PHASE EUNCTION WO ( I + XCOST f A ) 

iNClUEi'iF tihArt HAS THE UlKlCTluN OOSINL XIO 

JKAS6KHAK S MtThGu IS UStf), bUT rtlTH A COMPOSITE GAUSS gUAuRATURE 
JlSCRtTE UKUiiMATtS ARb STURLD I XM, WITH WEIGHTS IN A 

^ 4c X he zeroth ano first qkoek chahactekist ic kouts are stureu in XKO ANO XK 1 
4c#4T««T»)hE ZbROrH- ANO first UkOtR h-FUNCTlONS ARb SrURbO IN HMO AND HMl 
4 c4t4t*«*Th£ total reflectance k IS ALSO CuMPuTEU 

OlMENiluN A( 30)fXM( 3UI,XKO( 30),XRl( 30 I ,HMO( 30) t HMl ( 30 ) 

Pl = 3* l^li>yZ6536 
N-30 

i^4T»«4t.^(;ALCULATE THE N DISCRETE ORDINATES XMt AND WEIGHTS A 
call lNTRVL.(Nf XMfA) 

4c<c4c4t«*caMPUTb FOk i/AklOuS VALUES OF ANISOTROPY PARAMETER X, ANJ ALrTFOO WU 
00 20 0=1,3 $ X=2.0-1.0«J 

DO 20 K=l,2 $ WU= 1 .20-;20<^K 

4c*4t4c4c4tCALCOLATE THt ZEROTH ORDER H-FUNCTlON HMO, ANO THE CHARAC T E R I S T I C ROOTS 
XKO AT THE. DISCRETE UKulNATES XM 
L=0 

CALL HFUNC ( wQ , X , L , N , XM , A , X kO , HMO ) 

4c4t««**COMPOTE THE ZEROTH AND FlkST-GRUER MOMENTS OF H 
ALPHO=0. 

ALPH1=0. 

00 5 1=1, N 

ALPHO = ALPHOl-HMO( D^Ai I ) 

ALPH1 = ALPHI <-HMO( I )*A( 1 ) ♦XMC I ) 

5 CONTINUE 

C=X*WG*( U-WU) ^'ALPHl/IZ .-WU^ALPHO) 

PRINT 151 ,N,X, WO 

151 F0KMAT(*1 N=’«‘13« X=«F5.2* wO = *F5.2) 

PRINT 33 

33 F0RMAT(*0 MU HO(MU) Hl(MU)*) 

4c*«***CALCULATE THE FIRST-ORDER H-FUNCTION HMi AND CHARACTERISTIC ROUTS 
XKl AT DISCRETE ORDINATES XM 
L = 1 

CALL HFUNC ( WU , X , L , N , XM , A , XK 1 , HMI ) 

♦♦♦♦♦♦CALCULATE H-FUNCTIONS AT ARBITRARY Xl USING CHARACTERISTIC ROOTS 
DO 30 M=l,21 $ Xl = .05<^(M-l ) 

HO = I. 

Hl=L. 

00 35 1=1, N 
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APPENDIX A 


HU=^H0*(X1«-XM( I n/XM< n/( 1.>XK0( 1 J*Al ) 

Hl = Hl*(Xl«-XM(in/XM(l)/(l.«-XXl(II*Xll 

35 CONriNOE 

PRINT 34,Xi,M0fhl 
3^ fORMAT(F5.i:,tFi5.d) 

30 CUNTlNdt 

PRINT 36,ALPrl0tALPHl,C 

36 FOkMATI^O AL PHAO =*F 10 .6* ALPH A 1= oF 10 . 8« C = *Fl0.d) 

♦*****C0MPUTE REFutCTIuN FUR VARIUUS VALUES UF OIKECTICN JF INCIDENT LIGHT, XIO 
OU 13 1X0=1, 3 1 X10= l.i!-. 2*1X0 

PRINT 15 

® 15 FORMAT 1*1 X A'il MUG I HI T " -PHl - 10 II 

t I I/MUO*» 

HOO=l. 

H10=l. 

00 11 1=1, N 

HOO =H00 *(X104-XMmi/XM(H/(l.»XK0(l)*X10I 

11 HIO =hlO *( XiO«-Xf1( I » i/XM( 11/ 1 l.^XKK I ) *XlOi 

OU 10 1X=1,10 i Xl=CUS(PI*( IX-lt/18.1 

THEI =10.*( lX-1 ) 

♦ ♦***#tOMPOTE MU-UcPtNUcNT PAkT OF K EFL tC T I ON, EO 
HO = 1 . 

HI =1. 

00 12 1=1, N 

H0=H0*1 X1+XM( 1 J 1/XMI I J/( l.»-XKO( 1 I*X1 ) 

12 Hl = ril*I Xl*XM( lll/XM(I)/(l.tXKl(H*XU 
E0=w0*X10/(Xl+X10l/‘t.*H0*H00*( l.-C*lXl*X10I-X*l l.-W0I*Xl*X10» , 

♦♦♦♦♦♦CUMPOTc PHl-OEPENOENT PART OF RtFLECTluN, El 
00 16 IPH = 1,5 * PhI = Pl*< IPh-ll/',. 

PHI0=PHI*180./PI 

EliX^SOKTI I 1.-X1*X1»*( I.-XIO + X 101 1 *h 1 *H 10*CUS { PH I I * WO*X 10/ ( X I* XI 01 

*/A. . * 

ET=E0*E1 
ETX0=ET/X10 

PRINT 21,X,Ru,X10,THET ,PriIO,EO,El,tT,ETXO 
21 FOKMAT(3F/.2,F0.1,F7.1,2F10.5,F10.6,F10.4) . 

IFIABSIED.LT. I.e-IO.ANO.IPH.EO.IJ GU To 19 
16 CUNTINOE 

19 CONI INOE 

10 CONTINUE 

♦♦♦♦♦♦calculate total reflectance R 

K = 0. 

00 17 1=1, N 

1 7 K=K*HMO( U*( ( 1-X*( l-rtO»*XM( I) *X1 01 / 1 X 10«-XM( 1 1 I -C T *XM I I ) *A I I I 
R=nU*H00*R/2. 
print 18, K 

18 FORMAT I *0 R = *F9.6I 

13 CONTINUE 

20 continue 

STOP » END 
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APPENDIX A 


SU6RJUT1NE HEUNC ( WU t X . L t N ,XM, A>XK« H ) 

THIS SUbROUTlNfc COMPOTES CHANOKASEKriARS H rUHCriOHS FOR THE CASE 
***** WOU«-XCOSIH£TA) . L IS THE OHOER OF THE H FUNCTION. 

DIMENSION XMUt, All). XKIl), Hilt 
***** IF x=0 wE HAVE THE ISOTROPIC CASE. 

IF(X.EO.O..ANU.E.E(j.U GO TO 30 
***«**SN IS THE SIGN OF S ON THE LEFT SIDE OF EVERY POLE. 

SN= (X/AbSIXn**L 
******FIN0 THE N ROOTS. 

DO 19 J=1,N 

****** ihe LAST ROOT IS TREATED AS A SPECIAL CASE 
IFIO.EO.Nt GO TO 13 

******ThE ROOT IS LOCATED TO THE LEFT OF 1/XMIJ) AND TO THE RIGHT OF 1/XMU*1) 
X1=1./XM( J) 

DT=U1-1./XM( J*m/2. 

CALL ROOTIXl.DT.SNtNfXMf A.MOtLtX) 

♦*****stoke the rout in XR 

XKCJJ=X1 
GO TO 19 

13 CONTINUE 

******NUH FIND THE LAST ROOT. FIRST SET THE ROOT EQUAL TO ZERO AND CHANGE 

* LATER IF NECESSARY 
ARIJI^O. 

* *****COMP'OTE SO. THE VALUE OF S AT K=0. 

SO=-l. 

DO 14 1=1, N 
XMS=XM( 1I*XM( I I 

PSl=HO/2./(L*U*(l-L + ( 1-Z*U*X*IXMS*( l.-WO*( 1-LI I-LJ I 

14 SO=SO*A( I »*2.*PSI 

******!(: so IS EQUAL TO ZERO OR GREATER THAN ZERO, THEN THE ROOT LIES VERY 

* NEAR K=0. , SO LET IT STAY ZERO 

IFISD.GE.O.) GO TO 19 

** IF NOW SO IS LESS THAN ZERO WHEN SN IS -1, THE ROOT MUST 

* lie to the right of 1/XM(H 

IF(SN*SO.GT.O. 1 GO TO 15 ■ 

*♦ IF SO IS LESS than zero WHEN SN IS +1, THE ROOT MUST 

* LIE TO THE LtFT OF 1/XMIN) 

X1=1./XM(N1 

OT=Xl/2. 

call R0UT(X1,0T,SN,N,XM,A,WC,L,XI 

XMJI=X1 

GO TO 19 

15 X1=1./XM(1I 

0T=-Xl/2. , , 

CALL ROOnxl,OT,-SN,N,XM,A,WU,L,XI 
XK(J»=X1 

19 CONTINUE 

******COMPUTE THE H-FUNCTION 

DO 20 M=1,N i X1=XM(M) 

HtM)=l. 

DO 25 1=1, N 

H(M) = H(M)*(X1*XM( m/XMIIl/I l.*XK(I)*Xll 
25 CONTINUE 

20 continue 

RETURN 

30 CONTINUE 

******StT H=1 FOR THE ISOTROPIC CASE WHEN L=1 
DO 35 1=1, N 
XK(II=1./XH(1) 

HI I 1 = 1.0 
35 CONTINUE 

RETURN 
END 
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SUi3RUUTlNe INTKVLtNfXM, A) 

SUbRQUTlNc CUMPUTbS T Hb N UlSCKiiTb VALUES UF AM> AiMU WEIGHTS A 3Y 
^ BKbAKli'^G UP THE INTERVAL BtTwEEN 0 ANl> l» INTO N/2 SUB- 1 N T cK VAL S 

<= ANJ USING A 2-POlNT GAUSSIAN UlVlSluN H^ITHIN EACH SUB-INTEKVAL 

UlHEN^IUN XN( n ,A( 1) ,Xl 50) 
b = 2.0 
f1=N/2 
N = 2«.'1 

SaK3 = SUKU3.) 

♦♦♦♦♦^calculate Xf THE BOUNUAKiES UF THE SUBINTEKVALS 
UU 10 J=lfM 

10 X( J) = (FLUAT(J)/FLQAT(M) 

AL CULAT £ THE UlSCRtTE VALUES XM» ANU wElbHTS A 
Du 20 J=2fM 

A(2*J)=A(2*J-l )=(X(Ji-X(J-l) )/2. 

XM(2<^J-1) = (X(J)<’XU-1) )/2*-A(2*J)/Swk 3 
2 0 XM( 2« J ) = l XU) i-X( J-l ) )/2.^AI2*J )/SUK3 

A(iJ=A(2)=X(i)/2. 

XM( 1 ) = X( I ) /2.-AI 1) /SgR3 
X.'U2)=Xll)/2.<'A( 1 J/Swk3 

return 

END 


SUBROUTINE RUG T ( X 1 , DT , SN • N ,XM , A# WU» L , X J 
DIHENSIUN XMU ) , A< 1) 

♦♦♦♦♦♦THIS SUBROUTINE FINDS THE ROOT OF TH S-FUNCJiuN WHICH IS KNLWN TO 
♦♦ Lib BbTwEbN Xl ANU X1-2*UT 

♦♦♦♦♦♦THE SIGN OF SN IS twUAc TU THE SIGN OF S NE Ak Xl 
SU=SN 
1T = 0 

10 xi = xi-or 

1 T=IT4-1 

♦ ♦♦ INITIALIZE S AND T HE DbK I VAT I Vb SS 

S=-l. 

ss=o. 

♦♦♦ COMPUTE THE S-FUNCTION ANU THE DERIVATIVE OF THE S-FINCTION, SS 
DO 11 1=1 fN 
XMS = XM( I ) ♦XMi I ) 

D=l.-XMS^Xi^Xl 

♦♦♦ THE FUNCTION PSI DEPENDS ON THE QkDtK L 

PSi = wG/2./(L<-l ) ♦( l-L+( 1-2.«L) ♦X*(XMS^( I.-WO* ( l-L) )-L) ) 

Si=A( I )^2.^PS1/U 
SS=SS<-Sl/D^2.«XMS*Xl 

11 S=S+Sl 

♦♦♦♦♦♦HALVE THE INCREMENT 
DT=DT/2* 

S has changed sign then zero has been, crossed — CHANGE SIGN'UF INCREMENT 
IFIS^SO.LT.OU D1=-DT 
S0 = S 

fHb NEwIUn increment S/SS is SMALLER THAN DT, THEN THE NEW 

♦ INCREMENT IS THE NErtTuN INCREMENT 

IF(ABS(S/SS).LT.ABS(DI ) ) UT=S/SS 

♦♦ IF THE INCREMENT IS NOT SMALL COMPARED TU THE ROOT THEN RE-ITERATE 
IF(ABS(DT/XIJ.GT.1.E-13) GO -TO* 10 

♦♦ THE ROOT IS NOW RETURNED IN Xl 

Xl=Xl’DT 
KtTURN 

END . ■ ’ ' 
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APPENDIX B 


FLOW DIAGRAMS OF COMPUTER PROGRAMS 

This appendix -shows flow charts of the main computer program (program CHSKR) 
and subprc^rams (subroutine HFUNC and subroutine ROOT) given in appendix A. 


Flow Diagram of Program CHSKR 
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APPENDIX C 


DATA OUTPUT FOR COMPUTER PROGRAM 
This appendix shows the results of the prc^ram for various parameters. 


N= 30 X= 1.00 wU= 

. 

c 

c 

MJ 

HOC MO » 

HKMU) 

0.00 

i.oooooboo 

1.00000000 

.05 

1.13657325 

1.03582135 

• ip 

1 . 24/3 So 2 

1.05609911 

.15 

1. 35063305 

1.07122049 

.20 

i • 45U35 104 

1.06330765 

• 23 

1.54/32597 

1.09333557 

.30 

1.64252208 

1. 10165853 

.35 

i. 73640360 

1. 10922922 

.40 

1.82927553 

l.il5666ol 

♦ 

1.92134957 

1.12140966 

• 30 

2.012 //8/6 

1.. 12652112 

.53 

2.10367745 

1.13112161 

.60 

2.19413309 

1. 135268 3 7 

.65 

2.26421413 

1. 13908303 

. fO 

2.37397503 

1. 14255555 

.75 

2.‘t634596l 

i . 145 7469 3 

.60 

2. 552 7044 7 

1. 14669126 


2 . 641 73664 

1. 15141713 

.90 

2. 73056765 

1. 153946 73 

.9 5 

2.819272^2 

1. 15tj306 /d 

1.00 

2.90761073 

1. 1565C680 


AlPHAO=<J.OOOOOOOO AlPHA 1= 1. t = O.OOOOOOOU 


X 

wo 

MUO 

THtT 

PHI 

1.00 

1.00 

1.00 

0.0 

0.0 

1.00 

1.00 

1.00 

10.0 

0.0 

1.00 

l.OO 

1.00 

20.0 

0.0 

1.00 

1.00 

1.00 

30.0 

0.0 

1.00 

1.00 

1.00 

40.0 

0.0 

1.00 

1.00 

1.00 

50.0 

0.0 

1.00 

1.00 

1.00 

60.0 

0.0 

1.00 

1.00 

1.00 

70.0 

0.0 

1.00 

1.00 

1.00 

80.0 

0.0 

1.00 

1.00 

1.00 

90.0 

0.0 


R= 1. 000000 


10 

I 1 

1 

I /MUU 

1 .05.6^1; 

0.00000 

1.056920 

1.0569 

i.O:>616 

0. UOOOO 

1.055163 

1.0552 

1.0A<i/5) 

0.00000 

1.049752 

1 . 0498 

1.0A025 

0. 00000 

1.040251 

1.0403 


0.00000 

1.025826 

1.0258 

l.OObOS 

0.00000 

1.005051 

1.0051 

.97546 

0.00000 

.975463 

. .9755 

.9325)1 

0. 00000 

.932506 

.9325 

. U6609 

0. 00000 

. 666092 

.8661 

.72695 

0.00000 

.726953 

. 7270 
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APPENDIX C 


X 


MUO 

THLT 

PHI 

10 

1 1 

I 

I/hUO 

i .00 

1 .00 

.ao 

0.0 

0.0 


o.oaouo 

. 624753 

1.0309 

1.00 

1.00 

.60 

10.0 

0.0 


.01553 

.839611 

1.0495 

1 .00 

1.00 

.60 

10.0 

4 5.0 


.01098 

. 835063 

1.0436 

1 . 00 

1.00 

.60 

10.0 

90.0 

.02408 

-.00000 

.624083 

1.0301 

1.00 

1.00 

.60 

10.0 

135.0 

.82408 

-.01096 

.613103 

1.0164 

1.00 

1.00 

. 60 

10.0 

180.0 

. o 2408 

-.01553 

.606555 

1.0107 

1.00 

1 .00 

.80 

20.0 

0.0 

.82200 

.03132 

.653322 

1.0667 

1.00 

1.00 

.60 

20.0 

45.0 

.82200 

.02215 

. 844148 

1.0552 

1 . Oo 

1.00 

. 60 

20.0 

90.0 

.82200 

-.00000 

.822000 

1.0275 

1 .00 

1.00 

.80 

20.0 

135.0 

.82200 

-.02215 

.799851 

.9998 

1.00 

1 .00 

.60 

20.0 

160.0 

.82200 

-.03132 

. 790677 

.9885 

1 .00 

1.00 

• 60 

30.0 

0.0 

.81827 

.04767 

. 665940 

1.082 V 

1.00 

1.00 

. 60 

30.0 

45.0 

.81827 

.03371 

•651979 

1.0650 

1.00 

1.00 

• 60 

30.0 

90.0 

.81827 

-.00000 

•816273 

1.0226 

1.00 

1.00 

• 60 

30.0 

135.0 

.81827 

-.03371 

.764566 

. 980 / 

l.OO 

i.OO 

.60 

30.0 

180.0 

.81827 

-.04767 

. 770606 

.9633 

1.00 

1.00 

• 60 

40.0 

0.0 

. 812^5 

.06466 

.677327 

1 . 096 / 

1.00 

1.00 

.60 

40.0 

45.0 

.61245 

.04586 

.658325 

1.0729 

1.00 

1 .00 

.60 

40.0 

90.0 

.81245 

-.00000 

.812448 

1 . 015 O 

1.00 

1.00 

. 60 

40.0 

135.0 

.81245 

-.04586 

. 766572 

.9562 

1.00 

1.00 

.60 

40.0 

160.0 

.81245 

-.06468 

.747569 

.9345 

1.00 

1.00 

.60 

30.0 

0.0 

. 8 CJ 10 

.08333 

.687023 

1.1086 

1.00 

1.00 

. 60 

30.0 

45.0 

.80370 

.05892 

. 862617 

1.0783 

i.OO 

i.OO 

• Hi) 

50.0 

90.0 

.80370 

-.00000 

.603695 

1.0046 

1.00 

1.00 

.60 

50.0 

135.0 

.80370 

-.05892 

. 744773 

.9310 

1.00 

i.OO 

• 60 

50.0 

180.0 

.30370 

-.08333 

.720367 

.9005 

1.00 

1.00 

• 60 

bO.O 

0.0 

. 740^*7 

. 10345 

.893911 

1.1174 

1.00 

i.OO 

.60 

60.0 

45.0 

. 7904 7 

.07315 

• 66 3 6 1 3 

1.0795 

1.00 

1 .00 

• 60 

60.0 

90.0 

. 79047 

-.00000 

. 790466 

.9881 

1.00 

1.00 

. 60 

60.0 

135.0 

.79047 

. -.07315 

.717319 

. 8966 

1.00 

i.OO 

• 60 

60.0 

IdU.O 

. 79047 

-.10345 

.687021 

.8588 

1.00 

1.00 

• 60 

/ O.O 

0.0 

. 76 96 0 

.12568 

.895280 

1.1191 

1.00 

i.OO 

• 60 

70.0 

45 . 0 

. 76960 

.08887 

.658468 

1.0731 

1.00 

1.00 

.60 

70.0 

90.0 

. 76960 

-.00000 

. 769595 

.9620 

1.00 

1.00 

.60 

70.0 

135.0 

. 76960 

-.06687 

.660723 

.8509 

1.00 

1.00 

. 60 

70.0 

160.0 

. 76960 

-.12566 

.643910 

.8049 

1.00 

1.00 

• 60 

bO.O 

0*0 

. 73320 

.15019 

. 663396 

1.1042 

1.00 

1.00 

.60 

60.0 

45.0 

.73320 

•10620 

.639405 

1.0493 

1.00 

1.00 

.60 

60.0 

90.0 

. 73320 

-.00000 

.733204 

.9165 

1.00 

1.00 

• 80 

60.0 

135.0 

. 73320 

-.10620 

.627002 

.7638 

1.00 

1.00 

. 80 

60 . 0 

18 C .0 

. 73320 

-.15019 

.563011 

.7288 

1.00 

1.00 

.60 

90.0 

0.0 

.63816 

.17230 

.610460 

1.0131 

1.00 

1.00 

. 60 

90.0 

45.0 

.63818 

.12164 

.760013 

.9500 

1.00 

1.00 

• 60 

90.0 

90.0 

.63818 

-.00000 

.036176 

. 7977 

1.00 

1.00 

• 60 

90.0 

135.0 

.63818 

-.12164 

.516339 

• 6454 

1.00 

1.00 

• 60 

90.0 

180.0 

.63818 

-.17230 

.465872 

.5823 


R= 1. 000000 
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X 

wu 

rtUO 

THtl 

PHI 

10 

1 1 

I 

■ I/MUO 

I .00 

1.00 

.60 

0.0 

0.0 

.59814 

0. OOOOO 

. 598137 

.9969 

l.OO 

1.00 

• 60 

10.0 

0.0 

.59829 

.01728 

.615570 

1.0260 

l.OO 

1.00 

.60 

10.0 

45.0 

.59829 

.01222 

.610508 

1.0175 

1*00 

1.00 

.60 

10.0 

90.0 

.59629 

-.00000 

.598287 

.9971 

i.OO 

1.00 

.60 

10.0 

135.0 

.59829 

-.01222 

. 586065 

.9768 

1.00 

i.OO 

• 60 

10.0 

IdO.O 

.59829 

-.01728 

.581003 

.9083 

1.00 

1.00 

• 60 

^0.0 

0.0 

.59873 

.03498 

.633712 

1.0562 

1.00 

1.00 

• 60 

^0.0 

45.0 

.59873 

. Q 2473 

.623467 

1.0391 

1.00 

1.00 

• 60 

20.0 

90.0 

.59873 

-.00000 

. 598734 

.9979 

1.00 

1.00 

.60 

^0.0 

135.0 

.59873 

-.02473 

.574000 

.9567 

1.00 

1.00 

.60 

20.0 

180.0 

.59873 

-.03496 

.563755 

.9396 

1.00 

i.OO 

• 60 

iO.O 

0.0 

.59946 

.05354 

.653001 

1.0863 

i.OO 

1.00 

.60 

30.0 

95.0 

.59946 

.03786 

.637320 

1.0622 

1.00 

1.00 

• 60 

30,0 

90.0 

. 59946 

-.00000 

. 599463 

.9991 

1.00 

1.00 

.60 

30.0 

136.0 

.59946 

-.03786 

. 561606 

.9360 

1.00 

1.00 

• 60 

30.0 

100.0 

. 59946 

-.05354 

.545925 

.9099 

1.00 

1.00 

• 60 

^0.0 

0.0 

.60042 

.07351 

.673935 

1.1232 

1.00 

1.00 

.60 

40.0 

45.0 

.60042 

.05198 

.652404 

1.0873 

1.00 

1.00 

• 60 

' tO.O 

90.0 

.60042 

-.00000 

.600425 

1.0007 

1.00 

1.00 

.60 

40.0 

135.0 

.60042 

-.05198 

. 548445 

.9141 

i.OO 

1.00 

.60 

40.0 

180.0 

.60042 

-.07351 

.526915 

.8782 

1.00 

1.00 

.60 

50.0 

0.0 

.60148 

.09561 

.697088 

1.1618 

1.00 

1.00 

.60 

50.0 

45.0 

.60148 

.06761 

.669085 

1.1151 

1.00 

1.00 

• 60 

50.0 

9 C .0 

.60148 

-.00000 

.601479 

1.0025 

1.00 

i.OO 

.60 

50.0 

135.0 

.60148 

-.06761 

. 533873 

• 8898 

1.00 

1.00 

.60 

50.0 

180.0 

.60148 

-.09561 

.505870 

.8431 

1.00 

1.00 

.60 

60.0 

0.0 

.60222 

.12083 

. 723051 

1.2051 

i.OO 

1.00 

• 60 

oO . 0 

45.0 

.60222 

.08544 

.687661 

1.1461 

1.00 

1.00 

.60 

60.0 

90.0 

.60222 

-.00000 

.602223 

1.0037 

1.00 

1.00 

• 60 

60.0 

135.0 

.60222 

-.08544 

. 516786 

.8613 

i.OO 

1 .00 

• 60 

60.0 

180.0 

.60222 

-.12083 

. 481396 

.8023 

1.00 

1.00 

.60 

70.0 

0.0 

.60145 

.15059 

. 752041 

1.2534 

1.00 

1.00 

• 60 

70.0 

45.0 

.60145 

.10648 

. 707934 

1.1799 

1.00 

1.00 

.60 

70.0 

90.0 

.60145 

-.00000 

.601450 

1.0024 

1.00 

i.OO 

.60 

70.0 

135.0 

.60145 

-. 10 o 48 

.494966 

.8249 

1.00 

1.00 

.60 

70.0 

180.0 

.60145 

-.15059 

.450859 

. 7514 

1.00 

1.00 

.60 

60.0 

0.0 

. 59485 

. 18681 

• 781663 

1.3028 

1.00 

1.00 

• 60 

60.0 

45.0 

.59485 

.13210 

. 726946 

1.2116 

1.00^ 

1.00 

.60 

dO.O 

90.0 

.59485 

-.00000 

. 594849 

.9914 

1.00 

1.00 

.60 

bO.O 

135.0 

.59485 

-.13210 

.462 752 

.7713 

1.00 

1 .00 

.60 

ttO.O 

180.0 

. 59485 

-.18681 

.408035 

• 6801 

1.00 

1.00 

.60 

VO.O 

0.0 

.54853 

.22706 

. 775591 

1.2927 

1.00 

1.00 

.60 

90.0 

45.0 

.54853 

.16055 

. 709087 

1.1818 

1.00 

1.00 

.60 

90.0 

90.0 

.54853 

-.00000 

. 548533 

.9142 

1.00 

1.00 

.60 

90.0 

135.0 

.54853 

-.16055 

.387979 

.6466 

1.00 

1.00 

.60 

90.0 

180.0 

.54853 

-.22706 

.321476 

.5358 


K= 1.000000 
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N= 30 X= l.OO NU= .UO 


MU 

HUlMUl 

H 1 1 MU ) 

0.00 

1.00000000 

1.00000000 

.ot> 

1.0B746267 

1.02801646 

.10 

1.150122BA 

1.04362264 

.iu 

i.2037934d 

1.05515674 

.20 

1.25163311 

1.06431520 

.25 

1.29511719 

1.07187305 

.JO 

1.33511‘»07 

1.07826838 

.33 

1.37220207 
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.80 

20.0 

135.0 

.27658 

.01121 

. 287788 

.3597 

- 1.00 

.60 

.60 

20.0 

180.0 

.27658 

.01585 

.292430 

.3655 

- 1.00 

.60 

. 60 

30.0 

0.0 

.28095 

-.02423 

.256716 

.3209 

- 1.00 

.60 

• 80 

30.0 

45.0 

.28095 

-.01713 

.263814 

.3298 

- 1.00 

.60 

• 80 

30.0 

90.0 

.28095 

.00000 

.280949 

.3512 

-i .00 

.60 

• 80 

30.0 

135.0 

• 26095 

.01 713 

.296084 

.3726 

- 1.00 

.60 

.80 

60.0 

180.0 

.28095 

.02423 

.305181 

.3815 

- 1.00 

.60 

.80 

40.0 

0.0 

.28723 

-.03323 

.254003 

.3175 

- 1.00 

.60 

• 30 

hO.O 

45.0 

. 28 72 3 

-.02349 

.263735 

.3297 

- 1.00 

.60 

• 80 

40.0 

90.0 

.28723 

. 00000 

.287230 

.3590 

- 1.00 

.60 

.80 

40.0 

135.0 

.28 72 3 

.02349 

.310725 

.3884 

- 1.00 

.60 

• 80 

40.0 

180.0 

.28723 

.03323 

.320457 

.4006 

- 1.00 

.60 

.80 

50.0 

0.0 

.29554 

-.04315 

.252397 

.3155 

- 1.00 

.60 

. 80 

50.0 

45.0 

.29554 

-.03051 

.265034 

.3313 

- 1.00 

.60 

.80 

50.0 

90.0 

.29554 

.00000 

.295544 

.3694 

- 1.00 

.60 

• 80 

50.0 

135.0 

.29554 

.03051 

.326054 

.4076 

- 1 • Jo 

.60 

• 80 

50.0 

180.0 

.29554 

.04315 

. 338691 

.4234 

- 1.00 

.60 

.80 

60.0 

0.0 

. 30590 

-.05445 

.251448 

.3143 

- 1.00 

.60 

• 80 

60.0 

45.0 

.30590 

-.03850 

.267397 

.3342 

- 1.00 

. 60 

.80 

60.0 

90.0 

. 30590 

.00000 

.305901 

.3824 

- 1.00 

.60 

• 80 

60.0 

135.0 

.30590 

.03650 

• 344405 

.4305 

- 1.00 

.60 

.80 

60.0 

180.0 

.30590 

.05445 

.360354 

.4504 

- 1.00 

.'60 

.80 

70.0 

0.0 

.31 785 

-.06 769 

.249967 

.3125 

- 1.00 

. 60 

.80 

70.0 

45.0 

.31785 

-.04800 

.269850 

.3373 

- 1.00 

.60 

• 80 

70.0 

90.0 

.31785 

.00000 

.317852 

.3973 

- 1.00 

. 60 

• 60 

70.0 

135.0 

.31785 

.04800 

. 365854 

.4573 

- 1.00 

.60 

• 80 

70.0 

180.0 

•31 785 

.06789 

.385738 

.4822 

- 1.00 

.60 

.80 

80.0 

0.0 

.32913 

-.08485 

.244281 

.3054 

- l . OO ' . 

.60 ‘ 

. .80 

80.0 

45.0 

.32 913 

-.05999 

.269131 

.3364 

- 1.00 

.dO 

• 80 

80.0 

90.0 

. 32913 

.00000 

.329126 

.4114 

- 1.00 

.80 

.80 

80.0 

135.0 

. 32913 

.05999 

.389121 

.4864 

- 1.00 

.80 

.80 

80.0 

180.0 

. 32913 

.08485 

.413971 

.5175 

- 1.00 

.60 

.80 

' iO . 0 

0.0 

.32219 

-.11014 

.212043 

.2651 

- 1.00 

.80 

• 80 

*> 0.0 

45.0 

.32219 

-.07788 

.244303 

.3054 

- 1.00 

.80 

• 80 

90.0 

90.0 

. 32219 

.00000 

.322186 

.4027 

- 1.00 

.80 

.80 

90.0 

135.0 

.32219 

« 07 7 o 8 

.400069 

.5001 

- 1.00 

.80 

.60 

90.0 

180.0 

.32219 

.11014 

.432330 

.5404 


R= . 





APPENDIX C 


X 

^0 

MUO 

THET 

PHI 

10 

11 

1 

I/MUO 

- 1.00 

• dO 

.oO 

0.0 

0.0 

.21035 

0.00000 

.210346 

.3506 

- 1.00 

.dO 

.60 

10. 0 

0.0 

.21129 

-.00868 

.202406 

.3373 

-1 .uo 

.dO 

.60 

10.0 

46.0 

.21129 

-.00628 

.205009 

.3417 

- 1.00 

.dO 

• 60 

10.0 

90.0 

.21129 

.00000 

.211286 

.3521 

- 1.00 

.dO 

.60 

lO.O 

135.0 

.21129 

.00626 

.217568 

.3626 

- 1.00 

.dO 

• 60 

10.0 

IdO.O 

.21129 

•00686 

.220169 

.3669 

- 1.00 

.dO 

.60 

20.0 

0.0 

.21417 

-.01802 

. 196149 

.3269 

- 1.00 

.dO 

.60 

20.0 

45.0 

.21417 

-.01274 

.201427 

.3357 

- 1 .00 

.dO 

• 60 

20 . 0 . 

90.0 

.2141 7 

• 00000 

.214169 

.3569 

- 1.00 

.60 

• 60 

20.0 

135.0 

.2141 7 

.01274 

.225911 

.3782 

- 1.00 

.60 

.60 

20.0 

IdO.O 

.21417 

.01602 

.232189 

.3870 

- 1.00 

.dO 

.60 

30.0 

0.0 

.21915 

-.02771 

.191436 

.3191 

- 1.00 

.60 

.60 

30.0 

43.0 

.21915 

-.01960 

. 199553 

.3326 

- 1.00 

.dO 

• 60 

30.0 

90.0 

.21915 

•00000 

.219149 

.3652 

- 1.00 

.60 

.60 

3 0.0 

135.0 

.21915 

.01960 

.236746 

.3979 

- 1.00 

.dO 

.60 

30.0 

160.0 

.21915 

.02771 

. 246863 

.4114 

- 1.00 

. dO 

.60 

40.0 

0.0 

.22651 

-.03633 

•188180 

.3136 

- 1.00 

.dO 

.60 

40.0 

43.0 

. 2265 1 

-.02711 

.199408 

.3323 

- 1.00 

.dO 

• 60 

40.0 

90.0 

.22651 

.00000 

.226514 

.3775 

-l.Oo 

.60 

.60 

40.0 

135.0 

.22651 

.02711 

.253621 

.4227 

- 1.00 

.60 

• 60 

40.0 

IdO.O 

.22651 

.03633 

.264848 

.4414 

- 1.00 

.60 

.60 

60.0 

0.0 

.23669 

-. 0,5041 

• 186283 

.3105 

- 1.00 

.dO 

• 60 

60.0 

45.0 

.23669 

-.03565 

.201048 

.3351 

-1 . OO 

.60 

.60 

’ pO.O 

90.0 

.23669 

.00000 

. 236693 . 

.3945 

- 1.00 

.60 ^ 

.60 

60.0 

135.0 

.23669 

.03565 

.272338 

.4539 

- 1.00 

. 60 

.60 

60.0 

IdO.O 

.23669 

.05041 

.287103 

• 4785 

-l.OO 

.80 

.60 

6 0.0 

0.0 

.25026 

- • 064 76 

.185499 

.3092 

- 1.00 

.60 

.60 

60.0 

45.0 

.25026 

-. 045/9 

.204468 

.3406 

- 1.00 

.60 

• 60 

60.0 

90.0 

• 25026 

.00000 

.250263 

.4171 

- 1.00 

.60 

• 60 

60.0 

135.0 

.25026 

.04579 

.296058 

.4934 

- 1.00 

• dO 

.60 

60.0 

IdO.O 

.25026 

.06476 

.315027 

.5250 

- 1.00 

.60 

.60 

70.0 

0.0 

.26784 

-.08282 

. 185022 

.3064 

-l.OO 

.60 

• 60 

70. 0 

45.0 

.26784 

-.05856 

.209280 

.3486 

- 1.00 

.60 

• 60 

70.0 

90.0 

.26764 

.00000 

.267644 

.4464 

- 1.00 

• dO 

• 60 

70.0 

135.0 

.26764 . 

• 0585 o 

.326409 

.5440 

- 1.00 . 

.60 

.60 

. 70.0 

IdO.O 

.26784 

.08282 

.350667 

. 5844 . 

- 1 . 00 . . 

• dO. 

.60 

. tiO.O 

0.0 

.26927 

-.10746 

.181811 

.3030 

- 1.00 

.60 

.60 

60. 0 

45.0 

.20927 

-.07599 

.213285 

.3555 

- 1.00 

. 60 . , 

.60 

, 60.0 

90 . 0 

.26927 

.,00000 

.289270 

.4821 

- 1.00 

. 60 . 

• 60 

60. 0 

135.0 

.26927 

.07599 

.365256 ■ 

.6068 

- 1.00 

.60 

.60 

. 60.0 

IdO.O 

.26927 

.10746 

.396730 

.6612 

- 1.00 

.60 

.60 

VO.O 

0.0 

.30158 

-.14779 

. 153789 

.2563 

- 1.00 

.60 

.60 

90.0 

45.0 

.30158 

-.10451 

.197077 

.3285 

- 1.00 

.60 

.60 

90.0 

90.0 

.30158 

.00000 

.301582 

.5026 

- 1.00 

.60 

• 60 

90.0 

135.0 

.30158 

. 10451 

.406088 

.6 768 

- 1.00 

.60 

.60 

90.0 

IdO.O 

.30158 

.14779 

.449376 

.7490 


R = . 3952»9 


38 



REFERENCES 


1. Chandrasekhar, S,: Radiative Transfer. Dover Publ., Inc,, c.1960. 

2. Ralston, Anthony: A First Course in Numerical Analysis. McGraw-Hill Book Co., 

Inc., C.1965. 

3. Abramowitz, Milton; and Stegun, Irene A., eds.: Handbook of Mathematical Functions. 

Dover Publ., Inc., 1965. 


39 



TABLE I. - FUNCTIONS OBTAINED BY THE METHOD 

OF REFERENCE 1 AND PRESENT METHOD 




X = 

0 


X 

= 1 

M 

00 

o 

II 

o 

3 

Wo 

= 1.0 

Wo 

= 0.8 

Method of 
reference 1 

Present method 

Method of 
reference 1 

Present method 

Method of 
reference 1 

Present method 

0.05 

1.0820 

1.081 91 

1.1368 

1.136 57 

1.0876 

1.087 46 

.10 

1.1388 

1.138 81 

1.2474 

1.247 35 

1,1501 

1.150 12 

.20 

1.2286 

1.228 64 

1.4503 

1.450 35 

1.2516 

1.251 63 

.50 

1.4132 

1.413 27 

2.0128 

2.012 78 

1.4697 

1.469 71 

1.00 

1.5982 

1.598 22 

2.9078 

2.907 81 

1.7011 

1.701 11 


TABLED.- H^^) FUNCTIONS OBTAINED BY THE METHOD 
OF REFERENCE 1 AND PRESENT METHOD 



XWq = 1.0 

XO)o = 

0.8 

XCUo = - 

0.8 

X(Vq = - 

1.0 


Method of 
reference 1 

Present 

method 

Method of 
reference 1 

Present 

method 

Method of 
reference 1 

Present 

method 

Method of 
reference 1 

Present 

method 


1.0359 

1.035 82 

1.0281 

1.028 02 

0.9758 

0.975 82 

0.9702 

0.970 24 


1.0561 

1.056 10 

1.0436 

1.043 62 

.9637 

.963 71 

.9555 

.955 49 


1.0832 

1.083 31 

1.0642 

1.064 32 

.9488 

.948 78 

.9375 

.937 46 


1.1265 

1.126 52 

1.0966 

1.096 63 

.9277 

.927 66 

.9122 

.912 18 


1.1586 

1.158 52 

1.1201 

1.120 14 

.9137 

.913 72 

.8956 

.895 64 
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